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Abstract 

Over last decades, the study of laser fluctuations has shown that laser theory may be regarded as a 
prototypical example of a nonlinear nonequilibrium problem. The present paper discusses the fluctuation 
relations, recently derived in nonequilibrium statistical mechanics, in the context of the semiclassical laser 
theory. 

1 Introduction 

Nonequilibrium statistical mechanics aims at a statistical description of closed and open systems evolving 
under the action of time-dependent conservative forces or under time-independent or time dependent 
non-conservative ones. Fluctuation relations are robust identities involving the statistics of entropy 
production or performed work in such systems. They hold arbitrarily far from thermal equilibrium, 
reducing close to equilibrium to Green-Kubo or fluctuation-dissipation relations usually obtained in the 
scope of linear response theory |10l 1191 1241 1111 [HI E]. In a previous paper [2], we presented a unified 
approach to fiuctuation relations in classical nonequilibrium systems described by diffusion processes. We 
traced the origin of different fluctuation relations to the freedom of choice of the time inversion. The 
purpose of this paper is to illustrate the results of on the example of a phenomenological model of 
laser described by a stochastic differential equation. The semiclassical theory of laser describes the regime 
where, due to a large number of photons in the laser cavity, one may treat the electrical field classically, 
but the two level atoms are treated quantum mechanically [2UI [5] . The dynamical behavior of a single 
mode laser is then described by the equation of motion for the complex amplitude of the electric field Et: 

^ = {at- bEE)E, (1) 

where Et is the complex conjugate of Et- The function at is called the net gain coefficient and it takes 
into account the coherent emission and absorption of atoms and the losses. In the general case, at may 
have an explicit dependence on time, b is called the self-saturation coefficient. In most instances, it has a 
positive real part. There exist cases (with absorber) [21| where b has a negative real part, but we shall not 
consider them below. If the resonance frequency tJc of the laser cavity and the atomic frequency uja are 
exactly tuned then both at and b are real. In the case of detuning 22^, at and b are both complex. The 
equation of motion ([T]) describes the dynamical behavior of the laser field in a completely deterministic 
manner with the properties like coherence or spectral width lying outside the domain of the theory. The 
key to the understanding of such questions resides in the fluctuations of the electric field which are caused 
by random spontaneous atomic emissions. Such fluctuations may be accounted for by replacing Eq. ([T| 
by the stochastic differential equation 

^ = {at-bEE)E + r^{t,E,E), (2) 

with the noise r/^t, E, E) mimicking the effect of the random spontaneous emission of atoms in other 
modes, a purely quantum effect neglected in the semiclassical theory, but also the effect of vibrations of 
the cavity [20[ |9]. We shall take rjit, E, E) as a random Gaussian fleld with zero mean and delta-correlated 
in time. In the following, we shall look at two possible forms for 77, one additive and the other one 
multiplicative. The present paper is organised as follows. In Sect. 2, we recall the main results of [2]. 
In Sect. 3.1, we study the most elementary model of laser: the stationnary tuned laser with an additive 
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noise, and show that its dynamics satisfies the detailed balance. Sect. 3.2 is devoted to the fluctuation 
relations for a non-stationnary tuned laser. In Sect. 4.1, we examine the case of a stationnary laser with 
detuning. The detailed balance is broken here, but we show that its slight generalisation, the modified 
detailed balance, still holds. In Sect 4.2, we study the non-stationnary detuned case. In Sects. 5, we look 
at a slightly different case with the multiplicative noise. 

Acknowledgements. The author thanks Frangois Delduc and Krzysztof Gawedzki for encouragement, 
and Patrick Loiseau for his help in the numerical computation of Sect. 5.2. 



2 Fluctuation relation in diffusive systems [2] 

In [2], we dealt with arbitrary diffusion processes in J?"* defined by stochastic differential equation (SDE) 

X = ut{x) + Vt{x), (3) 

where x = ^ and, on the right hand side, Ut{x) is a time-dependent deterministic vector field (a drift), 
and vt{x) is a Gaussian random vector field with mean zero and covariance : 

(vlix)vi{y))=Sit~s)Dnx,y). (4) 

For the process solving the SDE ([sjl defined using the Stratonovich convention, we showed a detailed 
fluctuation relation (DFR) : : 

fio{dx)Po.T{x;dy,dW) e^-W) = /lUdy") Piriy' ; dx' , d{~W)), (5) 

where: 

• iJ,o{dx) = exp{—(poix))dx is the initial distribution of the original (forward) process, 

• jiQ^dx) = exp(— 93o(a;))da; is the initial distribution of the backward process obtained from the 
forward process by applying a time inversion (see below), 

• Po^t{x; dy, dW) is the joint probability distribution of the time T position Xj. of the forward process 
starting at time zero at x and of the functional [x] of the process (to be given later) that has the 
interpretation of the entropy production. 

• Pq^t{^\ dy, dW) is the similar joint probability distribution for the backward process. 

The key behind the DFR ([5| is the action of the time inversion on the forward system. First, the time 
inversion acts on time and space by an involutive transformation {t,x) — > {t* = T — t,x*). Second, to 
recover a variety of fluctuation relations discussed in the literature |14l 1151 HI [^ 1131 [T]. we allow for a 
non-trivial behaviour of the drift Ut under the time-inversion dividing it into two parts: 

ut — ut,+ + Ut,- (6) 

with Ut,+ transforming as a vector fleld under time inversion, i.e. Mfl +(a;*) — +{dkX*'''){x) Ut,^{x), and 
Ut,- transforming as a pseudo-vector field, i.e. itjl _(a;*) = —{dkx'''^){x) Ut,-{x). The random field Vt 
may be transformed with either rule: Vtl(x*) = ±{dkX*'^){x)vt (x). By definition, the backward process 
satisfies then the SDE 

X = u't{x) + Vt{x) (7) 

taken again with the Stratonovich convention. The functionnal which appears in the DFR depends 
explicitely on the functions ipo, (^o a-nd on the time inversion and has the explicit form: 

VKj, = A^^+ f Jtdt, (8) 
Jo 

where A^i/j = ipt{xt) — (po(a^o) with 

fij,{dx) = ejcp{—(pj,{x))dx = exp(— iy5o(3;*))cia;* = jj.'o{dx*), (9) 

and where 

Jt = 2ut,+ ■ dt^{xt){xt - ut,-{xt)) - V ■ ut,-{xt) (10) 

with dt(x) — Dt{x,x) and + = ul,+ — ^dyiDl^ {x,y)\y=x- The time integral in Eq. (jsj) should be taken 
in the Stratonovich sense. 
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The measures fiQ and /ig in the DFR ^ do not have to be normahzed or even normahzable. If they 
are, then distributing the initial points of the forward and the backward processes with probabiHties 
^o[dx) and ^'o(x), resperctively, we may define the averages 

{F)= J fio{dx) F[x], (F)' = J fi'oidx) E', F[x], (11) 

where Ex (-Ei) stands for the expectation value for the forward (backward) process atarting at x. From 
the DFR one may derive a generalisation of the celebrated Jarzynski equality [121 113j . 



(exp(-WV)) = l, (12) 

which may be viewed as an extension of the fluctuation-dissipation theorem to the situations arbitrarily 
far from the equilibrium. Note that the relation (12 1 implies the inequality (W^) > 0. 

To reformulate the DFR in a form where the entropic interpretation of is clearer, consider the 
probability measures M[dx] and M'[da::] on the spaces of trajectories of the forward and and of the 
backward process, respectively, such that 

(F) = J F[x]M[dx], {F}' = J F[x]M'[dx]. (13) 

The DFR may be reformulated in the Crooks form :5, as the identity 

{FeM~W^)) = {Fy, (14) 



where F[x] = F[x] with xt = a;T-t) s^nd the relation (14 1 implies the equality 

M'[dx] = exp(-Ty^[x]) M[dx], (15) 

for the trajectory measures with Al'[dx] — M'[dx]. By introducing the relative entropy S{M\M') — 
J ln( ^}^^^-^ ) A/ [da;] of the measure M' with respect to M, we infer that 

(W^) ^ S{M\M'). (16) 

Thus the inequality (Wj.) > follows also from the positivity of relative entropy. One may postulate that 
(Jt) dt describes the mean entropy production in the environment modeled by the stochastic noise: 

f {Jt)dt^A^Seuv. (17) 

Jo 

This is coherent with the previous result and particular cases, see [7l[8l[l6]. We may then interprete J^Jtdt 
as the fluctuating entropy production in the environment. An easy calculation leads to the relation 

(W^) =S(A^)-S(Mo) + A^Sen. +S'(ArlMr)> (18) 

where flt{dx) — ea;p(— i^t(x)) dx is the measure describing the time t distribution of the forward process if 
its initial distribution were /io(da;). S(p,t) = / (pt{x)fit{dx) is the mean instantenous entropy of the forward 
process Xt and S{fij.) — S{jj,o) is its change over time T. We could interprete ^pt{xt) as the fluctuating 
instantenous entropy. In general, p,^ is not linked to of formula (j9|. The relative entropy S{fir^\^j,) is 
a penalty due to the use at time T of a measure different than fi^. In the case where /t^ = (J'tj (W^t) 
the mean entropy production in the system and environment during time T and we could interpret Wj. 
as the corresponding fluctuating quantity. After a simple calculation [T^, one gets 



T rT 



A^Senv= {Jt)dt= dt [2ut^+{x)-dt\x){jt(x)dx-ut,^{x)fit{dx))-{V-ut,-){x)flt{dx)], (19) 
Jo Jo J 

S{fL^)-S{iio)^ [ dt [ jt{x)-V0t{x)dx, (20) 



where jt is the probability current at time t with the components 



ft = {ut - Id'/dj) exp(-(^t) (21) 



that satisfies the continuity equation 

dtexp{-(fit) + dift = 0. 
We shall apply now these results to three type of semiclassical single-mode laser. 
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3 Tuned laser with additive noise 



3.1 Stationnary case 

Let us consider the most common model of a stationnary laser with no detuning and with an additive 
form of the noise [20i i9]. Its dynamics is described by the SDE 

^ = (a - bEE)E + n, (22) 
with a and b real, 6 > 0, and with white noise rj with mean zero and covariance 

{■ntVt') = D5{t~t'), (23) 
iVtVt') = iVtVt') = 0. 

We can write the covariance matrix in the {E, E) space as 

d-D{\ I). (24) 



The equation ( 22 1 has then the form of the Langevin equation describing equilibrium dynamics of the 
process Et = (£t, St): 

^^-UvH^, + r^ (25) 

for Hab{E) = j^[b{EE)^ — 2aEE]. The Einstein relation is satisfied for the inverse temperature equal to 
1 implying that the Gibbs measure 

tJ.ab{dE) = Z-^^ exp{-Hab{E) dE (26) 

is invariant, has a vanishing probability current j, and satisfies the detailed balance 

Hab{dEo) Po.T{Eo;dE) = fiab{dE) Po.T{E;dEo). (27) 

This relation is a particular case of the detailed fluctuation relation ([5| where the time inversion acts 
trivially in the spatial sector, i.e. E* = E, the pseudo-vector part of the drift is taken zero, and we start 
with the Gibbs measure fiab for the forward and the backward processes. In this case both processes 
have the same distribution and = 0. The relation (271 may be projected to the one for the process 
It = EtEt describing the the intensity of the laser: 

f^abidlo) Po.TiIo;dI) = tlabidl) Po,T (I dig) . (28) 

The fluctuating entropy production in the environement may be identified with the heat production Q 
which is a state function here: 

A^g = / Jtdt ^ -H{E^) + H{Eo). (29) 
Jo 

This relations is the first principle of the thermodynamics in the case with no work applied to the system. If 
we start with the Gibbs density then the mean entropy production in the environment Aj.Senv = (A^Q) 
vanishes (191 as well as the instantaneous entropy production and Wj,. If the process starts with an 
arbitrary measure jj,o{dE) then at subsequent times the measure is 

At(dB) = J fioidEo) PoAEo;dE) (30) 

converging at long times to the invariant measure fj,ab{dE). During this process the mean rate of heat 
production (qt) in the environment is (19 1 



(qt) = (Jt) = - j (VHab ■ jt){E) dE . (31) 
After an integration by part, this may be written as 

{qt) = - J Hab{E) dtMdE) . (32) 
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3.2 Non-stationnary case 

3.2.1 Non-stationary net gain coefficient 

Let us consider now the SDE 

— = {at- bEE)E + r, (33) 
at 

with an exphcit time dependence for the (real) net gain coefficient at, with & > 0, and with the white 
noise rj as before. The explicit time dependance at may result from an external manipulation. In the 
matrix notation, the last equation takes the form 

riE 1 

-^-^dVHt + ,. (34) 

with Ht = Hatb- Here, we are outside the scope of the detailed balance and we enter in the world of 
transient fluctuation relations. To find an interesting DFR in this case, let us search for an appropriate 
time inversion. For example, we may impose that the backward process is still described by a Langevin 
equation but with the hamiltonian Ht{E) = Ht'{E*). By assuming a linear relation E* = ME and by 
transforming the drift with the vector rule, we obtain for the drift of the backward process the relation 

u't{x) = -^MdM^ (V //;)(£;). (35) 

To assure that MdM^ = d, we shaU take M = 1 or M = D'^d, i.e. E* = E or E* = E = {E,E). In 
these two cases, Ht{E*) = Ht{E) so that H'^{E) — Ht*{E) and the backward process satisfies the same 
SDE as the forward process but with the time-dependence of the Hamiltonian reparametrized. With this 
choices, a small calculation gives 



f Jtdt = - r VHt{Et)-dEt = -Ht{Et) + Hq{Eo)+ r{dtHt){Et 
J Jo Jo 



dt . (36) 



The first principle of thermodynamics implies then that {dtHt){Et) dt is the work performed on the 
laser during a time T. Starting from the Gibbs measure for the forward and the backward process, we 
obtain the relation 

'''^{dtat)Itdt, (37) 



Wt = -A^E + f {dtHt){Et)dt = -A^F - ^ f 
Jo ^ Jo 



where A^F = — Eg is the change of the Helmholz free energy Ft = — In J exp{—Ht{E))dE. The 
DFR ([sjl takes here the form 

fioidEo) Po,T{Eo;dE,dW) exp(-W) = fi^{dE) P;,^T{E*;dEo,d{-W)), (38) 

where /it denotes the Gibbs measure corresponding to Ht- In this case, there is a non vanishing entropy 
production in the environnement given by 

A^Senv = (A^Q) = f {Jt)dt= f dt f Ht{E){dtipt){E) oxp{-ipt{E))dE, (39) 
Jo Jo J 

where jit[dE) = exp{—ifi{E)) dE is the distibution of Et if Eo is distributed with the Gibbs measure 
^o[dE). Note that, in general, jit 7^ /it. The associated Jarzynski equality (121 takes the form 

T 

( exp [ - J{dtHt){Et) dt\ ) = exp(-AF), (40) 


that is, explicitly. 

In fact, there is an infinity of Jarzynski equalities that correspond to different splittings of the drift 



lit = — ^d'VHt into Ut,± parts. The peculiarity of the Jarzynski equality with the functionnal of (371 
is that upon its expansion to the second order in the small time variation at = a + ht with ht <^ a one 
obtains the standart ffuctuation dissipation theorem [101 1191 1241 llll [3] 

"''s 'h = 

for s < t, where ( ■ • ■ )q is the equilibrium average in the stationary state with h = 0. 
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3.2.2 External coherent field 

Another frequent way to induce a non-stationary behavior of the laser is to add an external coherent 
field at the laser frequency, modulated with a time-dependant amplitude which is injected into the 

cavity [10]. The gain and the self saturation of the laser depends now on the total field Et + E^^^ , but 
the losses depend just of Et, so the equation \22\ becomes: 



dE 



=l^a-b\E + Er'fj{E + E'")-aEr'+ri, (43) 

where a is the part of the dissipation in the net gain coefficient a. This equation takes for = Et + Ef^^ 
the form: 

= {a-b\E""f)E""-aE'''' + ^^+v. (44) 
Upon denoting —aE^^ H ^ — = /t, this may be rewritten as 

— = -^rfVH. + ,. (45) 

with 

In the case where Sf^' is not infinitesimal, we are outside the linear response regime, but the Jarzynski 
relation (|40[) is always true with 



dtHt{E'°') = --^ ((a/t)iJ'°* + {dtft)E"") . 
In the limit of infintesimal ft, this Jarzynski relation gives once again the fluctuation dissipation theorem 

m ■■ 

"^^^'^1 - ^d^iETAt)^, (46) 



Sfs U.o D 
5{At)\ _ 2 



Sfs \^^o D 

4 Detuned laser with additive noise 
4.1 Stationary case 

For the stationary case with no tuiung |22| . 

dE 



ds {ETAt)^ . (47) 



(a - bEE)E + ?7, (48) 



with a = a\+ ia2 and b = bi+ 162 complex, 62 > 0, and with covariance of the noise i] given by Eq. ( 23 1. 
The detuning destroys the Langevin form of the equation because the drift cannot be put any more in 
the form u — — ^Vi? but, instead. 
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U = -^VHa,b,+lDUVHa2b2, (49) 

with n = (J^i J). It is easy to see that the probability current of the Gibbs measure /laiti {dE) is 

jiE) = i Z-^\^UVHa^t2 exp{-Ha^b^) = Z-^\^ {-ib2E^E + ia2E, -ibiEE^ - iaiE) exp{-Ha^t^) (50) 

and that it is conserved: V ■ j =0 because H depends only on the intensity /. It follows that the measure 
/iai6i {dE) is preserved by the dynamics. We are in a steady state j6]. The detailed balance breaks down 
due to the non- vanishing of current j . It is replaced by the modified detailed balance: 

lia.bAdEo) Po.T{Eo;dE) = fia^b^dE) Po,T{E,dEo). (51) 

This relation, once again, implies a detailed balance for the process for intensity : 

Ataii,i(d/o) Po.T{Io;dI) = fia.bAdl) Po,T{I;dIo). (52) 



6 



The relation (51 1 is a particular case of the DFR ^ where the time inversion acts in the spatial sector 
as the complex conjugation E* = E, with the vector and pseudo-vector parts of the drift equal to 

u+ = {{ai - biEE)E, {ai -biEE)E) , U- = {i{a2 - b2EE)E, -i{a2 - b2EE)E). (53) 

Here again the backward process that we obtain with this choice of time inversion has the same distribution 
as the forward one and the heat production 



Jo 



(54) 



is a state function. If the forward and the backward processes are distributed initially with the Gibbs 
density exp{~Haibi) then, in average, there is no entropy production in environment 

A^5e„. = (A^Q) = fljt)dt = - [dt [vH^.b, ■ j{E)dE = - [dt f H^.b, ■ Vj{E)dE = (55) 
Jo Jo J Jo J 

and = 0. We have the usual features of equilibrium. 

4.2 Non-stationnary case 

Introduction of a time dependence of the net gain coefficient to the previous model leads to the SDE 

^ = (at - bEE)E + ri (56) 
with an explicit time dependence for the net gain coefficient at — ai,t +ia2,t and b = bi +ib2 with b2 > 0. 



Here, the fluctuation relation can be developed exactly as in Sect. 3. 2 but now (381 becomes for E* — E : 
fia,_obAdEo) Po,T{Eo\dE,dW) = Ha.^^b^idE) Po,T{E,dEo,d{-W)), (57) 
where fiai tbi denotes the Gibbs measure corresponding to Ha^ ^bi and 

^ -A^Fa.b, + f {dtHa,,b,){Et)dt ^ -A^Fa.b, - ^ [ {dtai,t)Itdt, (58) 
Jo ^ Jo 

The corresponding Jarzynski relation takes the form 



exp ll I {dta,,)Itdt]) = exp( "^'"^^^"^'° ) - 



compare to (41 1. The second order expansion in the small time variation at = a + ht with ht — /ii,t +ih2,t 



gives now the fluctuation-dissipation relations 

0/ll,s Ihso P* 0h2,s !ft=o 

see [3] for the details. 

5 Tuned laser with multiplicative noise 
5.1 Stationnary case 

It is not always clear a priori whether the noise is better represented by a multiplicative or additive model. 
In laser theory, when the randomness is due to pumping, it is more reasonable to use the multiplicative 
model of noise [20!. The stationary laser dynamics is then described by the non-Langevin SDE for the 
complex amplitude Et : 

^ = (a - bEE)E + ritE, (61) 

with a real, b positive and the white noise rit as before. In complex coordinates, the covariance matrix 
(|4| takes now the form 
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and, on the diagonal, 

d{E) = D{E,E) = DEE(^ ° J ). (63) 
One can show directly that the density exp(— (/:)(7)), where 

^(/) = ^/ + (l-^)ln/ (64) 

and / = EE, is preserved by the dynamics and corresponds to the vanishing current, leading to the 
detailed balance 

exp(-</j(£;„)) dEo Po,t{Eo; dE) = exp(-^(i;)) dE Po,t{E; dEo). (65) 

It is normalisable if a > 0. In this case, the normalized measure ii{dE) = exp{—tp{I)) dE is invariant 
and we are once again in an equilibrium case. There is no invariant probability measure when a < 0. 
Note the the intensity / satisfies here a closed SDE 

^ = 2(a - bl)l + (77t + 77t)J (66) 
that should be taken with the Stratonovich convention. 

5.2 Non-stationnary case 

Introduction of a time dependence of the net gain coefficient to the previous model results in the SDE 

^ = (at - bEE)E + r]tE. (67) 

With E* = E OT E* — E and the vector rule for the time-inversion of the drift, the backward process 
solves the same SDE with at and rjt replaced by at* and r]t* ■ This time reversal corresponds both to the 
so called reversed protocol and to the current reversal of the articles [l}|2]. The DFR ([5| takes now the 
form 

eM-MIo))dEo Po,T{Eo;dE,dW) exp(-H/) = exp(-vp^(/)) d£ P^,r(£*; d^, d(-W)) (68) 

with 

^*W = §/+(l-^)ln/ (69) 

and 

rT 



W^^ {dtipt){Et)dt = ~- {dtat)\nltdt. (70) 



The intensity process It satisfies the SDE (661 with the net gain coefficient a replaced by at- The 
backward intensity process is given by the same SDE with at and rjt replaced by at* and rjt' , leading to 
the DFR ([5} 

exp(-v5o(/o)))d/o Po,TiIo;dI,dW) exp(-M/) = exp(-<^T(/)) d/ Po,T{I\dIo,d{-W)). (71) 
Introducing the distribution of Wj, in the forward and the backward process by the relations: 

P,,iw)dw = I^M-Mio))dioPoMio-,di,dW)di 

J exp(-(y3o(/o)) dio 
P' mn,w - / exp(-yr(/o)) dh Pi,T{h; dl, dW) dl 

J exp(-v3T(io)) dIo 



we obtain by integration ( 68 1 the Crooks relation [3] 



Po,t{W) ^ Po^T{-W)exp{W - A^F) with A^F ^ Ft - Fo, (72) 

where Ft — — \nj exp{—ipt{I)) dl . In the case with positive ao and ar, we may derive the associated 
Jarzynski equality: 

(exp(-iyT)) = exp(-A^F), (73) 

where 

f exp(-u:>T(E))dE , , 

exp -A^F = \ ^\ ^ )'[ ,„ 74 

^ ^ ' Jexp{-<fo{E))dE ^ ' 
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or, explicitly 



cxp 



D 



(dtat) hi It dt ) = 



~2{aT-ao)/0 Y(2aT/D) 
d) r(2ao/D) ■ 



(75) 



Expanded to the second order in ht = at — a, the identity (731 induces the generalised fluctuation 
dissipation theorem (for a non-Langevin case): 



S {In It) \ 
6hs I 



— {In It ln/s)o 



(76) 



for s < t. Once again, it is the fluctuation dissipation theorem associated to the stochastic equation ( |67[ ), 
as it was demonstated in [3]. 



We did a numerical verification of the Crooks relation (721 for the case T — 1, at = 1 + t, b — 1 and 
D = 1. We realized with Patrick Loisearj^a Matlab computation. Below, we draw Po,i{W), Po.iiW) as 
a function of W and ln( "^"•^^'^^ ) as a function of — AiF. The simulation was done on 5000 initial 

conditions between and 10. For each initial condition, we considered 50 realizations of the noise. The 
interval of discretisation in time was 2~^^. 




Figure 1: Poj{W) as a function of W. Here (Wi) = 0.1023. 




Figure 2: j^{W) as a function of W. Here {W[) = 0.784. 



^Universite de Lyon, Ecolo Normalc Superieure de Lyon. 
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Figure 3: ln( 



Po.i(vy) 



) as a function of — Aii^. The continuum line is the identity function 



6 Conclusion 

We have discussed different fluctuation relations for a stochastic model of the semiclassical regime in 
a single mode laser. In particular, we showed that the stationary tuned laser with additive noise has 



an equilibrium state with detailed balance (27 1 and that the detuning preserves the features (51 1 and 



(551 of equilibrium. We also studied the non stationnary case, showing for the tuned and the detuned 



laser close to equilibrium the standart fluctuation-dissipation theorems ( 42 1 and ( 60 1 that extend to the 



appropriate Jarzynski equality (591 far from equilibrium. Finaly we studied laser with multiplicative 



noise. We specified in this case the detailed balance relation ( 65 I satisfied in the stationary case and 



the fluctuation-dissipation theorem (761. We also verified numerically the Crooks relation (721 in the 
non-stationary case. 
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